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Abstract 


Fall-winter  recurrence  of  current  reversal  from  westward  to  eastward  is  identified 
on  the  Texas-Louisiana  continental  shelf  using  the  current  meter  [Texas-Louisiana  Shelf 
Physical  Oceanography  Program  (LATEX-A)]  and  near  surface  drifting  buoy  [Surface 
Current  and  Lagrangian  Drift  Program  (SCULP-1)]  observations  in  1993  and  1994. 
Reversal  events  roughly  satisfy  the  Poisson  distribution  with  one  current  reversal  nearly 
every  twelve  days.  Synoptic  winds  seem  responsible  for  the  current  reversal  events.  Other 
processes  such  as  offshore  eddies  shed  from  the  Loop  Currents  and  river  run  off  are  less 
important  to  change  alongshore  flow  direction  at  synoptic  scales.  A  statistical  model  is 
established  to  predict  the  synoptic  current  reversal  using  the  surface  wind  observations. 
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1.  Introduction 


Bimodal  cycles  of  geopotential  anomaly  and  alternating  cyclonic  and  anti- 
cyclonic  regimes  are  important  dynamical  characteristics  of  the  Texas-Louisiana 
Continental  Shelf  (TLCS)  circulations  (Boicourt  et  al.  1998).  A  cyclonic  gyre  usually 
dominates  during  fall  and  winter  with  prevailing  westward  current  on  the  TLCS 
(Cochrane  and  Kelly  1986).  Satellite  infrared  images  and  ARGOS  drifters  in  spring  1989 
(Barron  and  Vastano  1994)  and  October  1992  (Vastano  et  al.  1995)  show  a  coherent,  ten- 
day  surface  current  response  over  the  entire  TLCS  to  the  establishment  of  a  high 
atmospheric  pressure  system  over  the  Gulf  of  Mexico.  Alongshore  westward  current  may 
stop  with  possibility  of  reversal  into  eastward  current  when  northerly  winds  are  sustained 
over  a  few  days  (Nowlin  et  al.  1998)  or  when  an  anti-cyclonic  circulation  regime  replaces 
the  dominating  cyclonic  regime  over  the  TLCS  (Oey  1995;  Li  et  al.  1997;  Cho  et  al. 
1998;  and  Ohhnann  et  al.  2001). 

Questions  arise:  Do  the  TLCS  current  reversals  (from  westward  to  eastward) 
occur  on  the  synoptic  timescale  that  dominates  the  surface  wind  forcing?  If  so,  what  are 
the  characteristics  of  the  TLCS  current  reversals?  What  is  the  major  forcing  to  cause  the 
synoptic  TLCS  current  reversal?  To  answer  these  questions,  establishment  of  temporally 
varying  current  velocity  data  with  high  spatial  resolution  is  needed. 

In  this  study,  the  optimal  spectral  decomposition  (OSD)  method  (Chu  et  al. 
2003a,  b)  is  used  to  process  the  current  velocity  data  collected  from  the  Texas-Louisiana 
Shelf  Physical  Oceanography  Program  (LATEX-A)  and  first  segment  of  the  Surface 
Current  and  Lagrangian  Drift  Program  (SCULP- 1)  into  temporally-spatially  gridded  data. 
With  this  data  set,  the  TLCS  current  structure  and  reversal  may  be  identified.  In 
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conjunction  with  the  local  surface  wind  data,  the  wind  effect  on  the  TLCS  synoptic 
current  reversal  may  also  be  identified. 

2.  Data 

Data  for  the  TLCS  fall/winter  current  reconstruction  were  collected  from  31  near¬ 
surface  current  meter  moorings  during  LATEX-A  at  10-15  m  depths  (Cho  et  al.  1998) 
from  April  1992  to  December  1994  (Fig.  1)  and  from  drifting  buoys  deployed  during 
SCULP-1  (Ohhnann  et  al.  2001)  from  October  1993  to  July  1994.  The  current  meter 
moorings  sample  the  current  velocity  (speed  and  direction),  temperature,  and  salinity 
every  5-min  to  2-hour  (mostly  30  min).  The  drifting  buoys  record  the  location  at  various 
times  with  inhomogeneous  area  coverage.  The  data  were  interpolated  into  a  uniform 
temporal  grid  with  time  step  of  3  hours. 

Meteorological  data  are  from  7  buoys  (represented  by  squares  in  Fig.  1)  from  the 
National  Data  Buoy  Center  (C-MAN)  in  the  LATEX-A  area.  Wang  et  al.  (1998)  have 
estimated  the  zero-crossing  spatial  scales  of  winds  (-350  km)  and  other  meteorological 
parameters  over  the  TLCS  shelf  after  analyzing  extensively  the  meteorological  data 
collected  during  the  LATEX-A  Program  observation  period.  Since  the  scale  of  the 
atmospheric  systems  is  larger  than  the  scale  of  the  continental  shelf,  the  horizontal  mean 
wind  speed  and  direction  from  7  buoys  are  used  for  analysis. 

3.  OSD  Method 

A  low-pass  (40  hr)  filter  is  used  to  get  rid  of  energetic  shelf  oscillations  caused  by 
tides,  inertial,  and  sea  breeze  from  the  original  current  mooring  and  drifter  observational 
data  (Chen  et  al.  1996;  DiMarco  et  al.  2000).  The  filtered  data  are  reconstructed 
(0.1°x0.1°  grid)  using  the  OSD  method  (Chu  et  al.  2003  a,  b). 
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Between  April  1992  and  December  1992,  the  LATEX-A  current  meter  data  have 
good  temporal-spatial  coverage,  but  the  SCULP- 1  drifter  data  are  not  available.  To 
reduce  the  reconstruction  error  caused  by  spatial  sparseness  and  temporal  interruption  in 
the  observational  data  (8  times  daily),  three  periods  (248  days)  are  selected  with 
reasonable  reconstructed  near  surface  circulation:  (a)  from  January  21  to  May  21,  1993 
(121  days,  968  observations),  (b)  from  December  19,  1993  to  February  17,  1994  (61 
days,  488  observations),  and  (c)  from  October  5  to  November  29,  1994  (56  days,  448 
observations). 

4.  OSD  Reconstruction  Skill 

Reconstruction  of  currents  from  sparse  and  noisy  observational  data  is  an  ill- 
posed  problem  that  may  not  have  a  unique  solution.  Theoretical  and  practical  evaluations 
are  conducted.  The  theoretical  evaluation  is  to  check  if  the  solution  is  robust  to  sample 
size  and  noise  level  and  asymptotically  stable,  i.e.,  if  the  reconstructed  solution 
homogeneously  converges  to  an  optimal  solution  when  the  number  of  observations  tends 
to  infinity  (Ivanov  et  al.  2001;  Chu  et  al.  2003a,  b). 

For  evaluation,  the  LATEX-A  data  before  the  reconstruction  are  divided  into 
testing  and  working  data.  The  testing  data  are  data  collected  at  four  current  meter 
stations:  2,  4,  18  and  24  (Fig.  1).  The  working  data  are  data  collected  at  all  the  other 
current  meter  stations.  The  working  data  are  used  for  OSD  reconstruction;  and  the  testing 
data  are  used  for  OSD  verification. 

Let  [u\rec\v\rec)]  and  [u\obs) ,v\"bs)\  be  i-th  reconstructed  and  observed  zonal  and 
latitudinal  velocity  components  (after  40  hr  low-pass  filter)  at  the  testing  stations.  Fig.  2 
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shows  [u\rec\v\rec)]  (dashed  curve)  and  [u\"bs) ,v\"bs)\  (solid  curve)  at  the  four  testing 
stations.  The  relative  root  mean  square  (RMS)  errors, 


mnse(w)  = 


rrmse(v)  = 


are  (17.1%,  7.3%)  at  Station-2,  (3.4%,  13  %)  at  Station-24,  and  are  larger  at  Stations-4 
and  Station- 18  (Table  1). 

5.  Fall-Winter  Synoptic  Current  Reversal 

Trajectories  of  77  SCULP-I  drifting  buoys  show  current  reversals  during  January 
2-7,  1994  (Fig.  3).  The  OSD  reconstructed  horizontal  velocity  vector  field  shows 
dynamical  characteristics  of  the  synoptic  current.  The  flow  pattern  on  30  December  1993 
is  characterized  as  a  cyclonic  gyre  with  a  strong  westward  TLCS  flow  as  the  northern 
flank  (Fig.  4a).  This  current  completely  reverses  to  the  east  on  3  January  1994  (total 
reversal,  Fig.  4b),  and  the  cyclone  re-occurs  (6  January  1994)  with  the  westward  current 
at  the  north-central  shelf  and  the  eastward  currents  at  the  western  shelf  and  shelf  break 
(partial  reversal,  Fig.  4c).  The  total  reversal  has  a  width  up  to  200-m  isobath  (depth  for 
shelf  break),  and  the  partial  reversal  has  a  width  smaller  than  200-m  isobath.  During  that 
period  (30  December  1993  to  6  January  1994),  southwesterly  winds  prevail,  and  no 
evident  offshore  eddy  appears  at  the  shelf  edge. 

Recurrence  of  the  TLCS  current  reversals  is  estimated  using  the  OSD 
reconstructed  velocity  fields  (230  days  of  duration).  A  period  T  (from  5  to  40  days)  is 
selected  as  the  time  window.  The  number  of  total  (or  partial)  reversals  is  counted  in  each 
T  day  window.  Let  (n0,  ni,  n2,  ...)  be  the  numbers  of  0-current  reversal,  1-current 


6 


reversal,  2-current  reversals,  and  m  be  the  all  realizations  for  a  given  T  day  window.  The 
probabilities  for  0-,  1-,  and  2-current  reversals  are  calculated  by 

P0(T)  =  ^,  PX{T)  =  ±,  P2 (T)  =  — ,  (1) 

in  m  m 

where  T’o(T),  T’i(T),  and  P2( T)  depend  on  the  period  T  (Fig.  5).  With  T  larger  than  20 
days,  the  probability  for  zero  current  reversal  is  less  than  0.2.  With  T  of  15  days,  the 
probability  for  one  reversal  reaches  0.5.  The  probability  PilT)  is  fitted  to  the  Poisson 
distribution  (Fig.  5) 

JPi(T)  =  -i(//r)A'exp(-//r),  *=0,1,2  (2) 

where  /u  =  0.08  day'1  is  the  mean  rate  of  current  reversal.  The  probability  of  time  interval 
between  successive  reversals  (r)  is  estimated  using  the  Poisson  distribution  (2) 
(Tikhonov  and  Khimneko  1998), 

p(T)  =  //exp(-//r) ,  (3) 

which  shows  that  the  synoptic  reversal  is  not  a  rare  event.  Its  characteristics  can  be 
obtained  from  analysis  of  temporal  and  spatial  current  variability. 

6.  Empirical  Orthogonal  Functions 
The  empirical  orthogonal  function  (EOF)  method  (Preisendorfer  1988)  is  used  to 
analyze  the  temporal  and  spatial  variability  of  the  three  sub-datasets  of  OSD 
reconstructed  velocity.  The  first  EOF  (Table  2)  is  able  to  account  for  80.2%  variance  for 
the  first  period  (January  21  -May  21,  1993),  77.1%  for  the  second  period  (December  19, 
1993-  February  12,  1994),  and  74.4%  for  the  third  period  (October  5  -  November  29, 
1994).  Therefore,  the  velocity  field  can  be  roughly  represented  by 

u(x,  y ,  t)  =  u(x,  y)  +  Ax  (t)u,(x,  y),  (4) 
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for  each  period.  Here,  u(x,  y)  is  the  mean  velocity  (background  circulation),  u,(x,v)  is 
the  EOF-1  mode  (nondimensional),  and  A\(t)  is  its  principal  component. 

The  background  velocity  field  u(x,y)  presents  a  cyclonic  circulation  pattern  (Fig. 
6)  with  mean  westward  TLCS  currents  similar  to  Cho  et  al.’s  (1998)  results.  During  the 
first  period,  cyclogenesis  over  the  shelf  produced  a  strong  winter  stonn  (then-called 
"Storm  of  the  Century")  of  12-16  March  1993,  which  strongly  affected  currents  during 
this  month.  The  mean  pattern  shown  in  Fig.  6a  strongly  resembles  the  cyclonic 
atmospheric  pattern  during  the  time  of  the  stonn. 

The  EOF-1  mode  presents  an  anticyclonic  circulation  pattern  for  positive  A\(t) 
(Fig.  7)  with  mean  eastward  TLCS  currents.  Comparison  is  conducted  between  the 
reconstructed  fields  from  two  independent  observational  data  sets  collected  during  the 
same  period.  High  coherence  (conelation  coefficient  ~  0.92)  is  found  in  A\{t)  from  19 
December  1993  to  12  February  1994  between  using  the  LATEX- A  current  meter 
mooring  data  (solid  curve)  and  the  SCULP-I  drifter  (dashed  curve)  data  (Fig.  8).  This 
shows  the  reconstruction  skill  of  the  OSD  method. 

The  TLCS  current  reversal  depends  on  the  relative  magnitudes  of  u(x,y)  and 

Aj(t) u,(x,y)  [positive  A\(t)\.  Since  u(x,y)  and  ut(x,y)  are  totally  opposite  on  the 

TLCS,  total  and  partial  current  reversals  can  be  simply  represented  by  two  thresholds  St 
and  Sp.  The  background  current  velocity  field  u(x,y)  changes  from  one  period  to  the 
other  (Fig.  6),  the  thresholds  for  the  reversal  also  vary.  For  the  first  period  (January  21  - 
May  21,  1993),  the  two  thresholds  are  given  by 

S.  =  23  m/s,  Sn  =  8  m/s  . 

i  p 
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Eight  total  reversal  events  are  identified  during  the  first  period  (Fig.  9a):  February  7-9; 
12-13,  March  5-1 1,  14-16;  March  31 -April  3,  April  9-12,  14-17,  and  May  9-12,  1993. 

Temporally  varying  weight  ratio  between  square  of  the  PC-1  to  sum  of  the 
squares  of  PC-2  to  PC-6  (Fig.  9b), 

7(0  =  [4(0]2/Zt4,(0]2,  (5) 

n= 2 

represents  importance  of  the  EOF-1  mode  relative  to  higher  EOF  modes.  The  total  and 
partial  current  reversal  events  correspond  to  large  77  values. 

7.  Wind  Driven  Mechanism 

Wind,  freshwater  influx,  and  offshore  eddies  affect  TECS  circulations  (Boicourt 
et  al.  1998).  What  is  the  major  forcing  for  the  TECS  current  reversal?  Analysis  of 
reversal  events  for  the  first  period  (December  1993  -May,  1993)  is  taken  as  a  sample  for 
the  study.  The  horizontally  averaged  surface  zonal  and  latitudinal  wind  components 
[(7(0,  V(t)\  (8  times  daily)  are  calculated  from  the  seven  meteorological  buoys.  The 
Texas-Louisiana  coastline  is  slanted  roughly  17°  to  the  zonal  direction.  The  alongshore 
wind  [(7ais(0]  is  the  projection  of  the  surface  wind  vector  to  the  axis  rotated  17°  anti¬ 
clockwise  from  the  latitude.  Comparing  Ua\s(t)  to  A\(t)  (Figs.  9a  and  9c),  the  total 
reversals  on  TLCS  are  associated  with  the  eastward  alongshore  winds  with  wind  speeds 
larger  than  2  m/s.  The  correlation  coefficient  between  (7ais(t)  and  A \(t)  is  0.82  for  the  first 
period  (sample  size:  968),  0.83  for  the  second  period  (sample  size:  488),  and  0.74  for  the 
third  period  (sample  size:  448).  For  significance  level  a  =  0.005,  the  critical  value  of  the 
sample  correlation  coefficient  is  0.254  for  sample  size  of  102,  and  is  less  than  0.254  for 
sample  size  larger  than  102. 
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Connection  between  t/ais(0  and  A  \(t)  (sampling  rate  of  8  times  per  day)  is  further 
investigated  using  their  temporally  varying  spectra  (Fig.  10)  from  the  Morlet  wavelet 
transform  (Morlet  et  al.  1982), 

O(t)  =  exTp(imt  - 12  /  2),  m  =  6,  (6) 

which  are  similar  at  frequencies  between  0.2  and  0.05  cpd.  This  may  confirm  the 
connection  between  synoptic  current  reversals  and  alongshore  surface  winds. 

8.  Statistical  Prediction  Model 

For  predicting  synoptic  current  reversal,  the  observed  surface  wind  [U(t),  V(t)] 
(not  alongshore  and  cross-shore  components)  is  used  with  a  simple  statistical  model, 

A  (0  =  a  [uit) -  u] + p  [V(t) -  v] + r ,  (7) 

where  U  and  V  are  temporally  averaged  wind  speeds;  a,P  and  y  are  regression 
coefficients.  The  regression  coefficients  are  calculated  from  the  data  collected  during  the 
second  period  (December  19,  1993  to  February  17,  1994), 

a  =  9.032,  P  =  3.754,  y  =  0.209  m/s.  (8) 

Nonzero  regression  coefficient  (y  *  0)  suggests  that  the  wind  forcing  may  not  be 
the  only  factor  causing  the  synoptic  TLCS  current  reversal.  This  is  because  the  current 
reversal  may  occur  ( At  ^0)  even  if  there  is  no  wind  anomaly  [U(t)  =  U,  V(t)  =  V], 
Since  the  two  regression  coefficients  for  the  wind  speeds  are  not  small  (a  =  9.032, 
P  =  3.754),  contribution  of  y  (-0.2  m/s)  to  A  \(t)  is  negligible  when  the  surface  winds 
are  stronger  than  0.4  m/s. 

The  Mississippi-Atchafalaya  River  discharge  increases  the  cross-shelf  salinity 
gradients.  With  salinity  increasing  offshore,  the  cyclonic  circulation  enhances  (Chen  et  al. 
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2000)  and  stronger  eastward  (alongshore)  winds  are  required  to  reverse  the  mean 
westward  TLCS,  i.e.,  the  threshold  of  reversal  is  higher. 

The  offshore  eddies  of  the  Loop  Current  origin  may  also  affect  the  synoptic 
current  reversal.  To  investigate  such  an  effect,  the  SCULP-I  drifter  data  are  reconstructed 
for  weak  alongshore  winds  (£4ig  around  0.4  m/s,  Fig.  11)  and  strong  alongshore  winds 
(C/aig  >  3.5  m/s,  Fig.  12).  For  weak  alongshore  winds,  offshore  eddies  drive  strong 
fluctuations  of  the  alongshore  flow  across  the  shelf,  and  interact  with  the  alongshore  jet 
(Fig.  1  lb).  For  strong  alongshore  winds,  the  offshore  eddies  appear.  They  can  only  cause 
a  local  acceleration  or  deceleration  of  alongshore  flow  but  cannot  change  its  direction 
[small  value  of  y  in  (7)]. 

For  strong  alongshore  winds,  offshore  eddies  do  not  appear  on  the  sea  surface 
(Fig.  12a)  but  occur  on  the  upper  slope  as  counter-rotating  pairs.  Such  offshore  eddies 
may  deflect  the  eastward  currents  into  southeastward  currents  in  the  western  TLCS  (Fig. 
12b). 

The  statistical  prediction  model  (7)-(8)  was  established  using  the  data  collected 
during  the  second  period  (December  19,  1993  to  February  17,  1994).  The  data  collected 
during  the  third  period  (October  5  to  November  29,  1994)  are  taken  as  the  independent 
data  set  to  verify  the  statistical  model.  Fig.  13  shows  the  comparison  between  the 
reconstructed  (from  the  third  period  data,  solid  curve)  and  predicted  (dashed  curve) 
principal  component,  A\{t).  The  correlation  coefficient  between  the  two  is  0.7.  The 
statistical  model  (7)  predicted  the  occurrence  of  the  three  total  reversal  events:  October 
14,  1994;  October  23,  1994  and  November,  1994  and  the  three  partial  reversal  events: 
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November,  2  1994;  November,  7  1994  and  November,  23  1994.  But  it  failed  to  predict 
two  partial  reversal  events:  November,  4  1994  and  November,  10  1994. 

9.  Conclusions 

(1)  The  OSD  method  is  used  to  analyze  fall-winter  LATEX-A  and  SCULP- 1 
data.  Theoretical  and  practical  evaluations  show  a  high  reconstruction  skill  of  the  OSD 
method.  The  mean  OSD  reconstructed  horizontal  velocity  field  shows  a  cyclonic 
circulation  pattern  near  the  Taxes-Louisiana  coast  with  the  westward  flowing  TLCS 
currents. 

(2)  Relative  RMS  errors  between  reconstructed  and  observed  (not  used  for  the 
reconstruction)  zonal  and  latitudinal  current  velocity  data  are  used  to  measure  the  OSD 
reconstruction  skill.  At  Station-24,  the  relative  RMS  error  is  3.4%  for  zonal  current  and 
13  %  for  latitudinal  current.  The  reconstructed  velocity  fields  are  consistent  using  the 
LATEX-A  current  meter  mooring  data  and  the  SCULP-I  drifter  data  collected  at  the  same 
period  with  high  coherence  (correlation  coefficient  ~  0.92)  4 ,(/)  from  19  December 
1993  to  12  February  1994.  This  also  shows  a  high  reconstruction  skill  of  the  OSD 
method. 

(3)  Temporally  varying  OSD  reconstructed  velocity  fields  show  that  the  synoptic 
TLCS  current  reversal  is  not  a  rare  event.  Probability  for  its  recurrence  satisfies  the 
Poisson  distribution.  With  period  larger  than  20  days,  the  probability  for  zero  current 
reversal  is  less  than  0.2.  With  period  of  15  days,  the  probability  for  one  reversal  reaches 
0.5. 

(4)  The  EOF-1  mode  (anticyclonic  circulation  pattern)  accounts  for  80.2% 
variance  for  the  first  period  (January  21  -May  21,  1993),  77.1%  for  the  second  period 
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(December  19,  1993-  February  12,  1994),  and  74.4%  for  the  third  period  (October  5  - 
November  29,  1994).  Occurrence  of  fall-winter  total  and  partial  synoptic  current 
reversals  on  the  TLCS  (i.e.,  altering  the  mean  cyclonic  circulation  pattern)  depends  on  the 
value  of  the  principal  component  corresponding  to  EOF-1.  When  the  value  exceeds  a 
threshold,  the  synoptic  current  reversal  occurs  on  TLCS.  Eight  total  reversal  events  are 
all  identified  during  the  first  period  by  this  method:  February  7-9;  12-13,  March  5-11,  14- 
lb;  March  31-April  3,  April  9-12,  14-17,  and  May  9-12,  1993. 

(5)  Alongshore  wind  forcing  is  the  major  factor  causing  the  synoptic  current 
reversal.  Other  factors,  such  as  the  Mississippi-Atchafalaya  River  discharge  and  offshore 
eddies  of  Loop  Current  origin,  may  affect  the  reversal  threshold,  but  can  not  cause  the 
synoptic  current  reversal. 

(6)  A  simple  statistical  model  is  developed  to  predict  the  total  and  partial  reversals 
of  the  TLCS  currents  using  the  surface  winds.  The  model  evaluation  using  independent 
data  set  shows  high  capability  to  predict  the  total  synoptic  reversals,  but  low  capability  to 
predict  the  partial  reversals. 

(7)  The  Texas-Louisiana  shelf  is  strongly  bimodal  having  a  summer  and  non¬ 
summer  seasonal  oscillation.  For  the  summer  season,  the  predominant  current  direction  is 
eastward,  which  is  opposite  to  the  winter  season.  Therefore,  only  the  fall-winter  data  are 
analyzed  in  this  study.  A  similar  analysis  may  be  used  to  see  if  synoptic  TLCS  current 
reversal  occurs  in  summer. 
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Table  1.  Relative  root  mean  square  errors  at  the  four  testing  stations 


Mooring 

Station 

rrmse(u) 

(%) 

rnnse(v) 

(%) 

2 

17.1 

7.3 

4 

14.9 

21.4 

18 

11.8 

27.6 

24 

3.4 

13.0 

Table  2.  Variance  of  the  first  six  EOFs. 


EOF 


01/21/93-05/21/93 


Variance  (%) 


12/19/93-04/17/94  10/05/94-11/29/94 


Fig.  1.  Geography  and  topography  of  Texas-Louisiana  continental  shelf,  LATEX-A 
current  meter  stations  (represented  by  the  symbol  ‘  •  ’),  and  meteorological  buoys 
(represented  by  the  symbol  ‘  □  ’). 
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Fig.  2a.  OSD  reconstructed  (dashed  line)  and  observed  (black  line,  after  40  hr  low-pass 
filter)  zonal  (upper  panel)  and  latitudinal  (lower  panel)  velocity  components  at  the 
stations  not  used  for  the  reconstruction:  Station-2. 
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Fig.  2b.  OSD  reconstructed  (dashed  line)  and  observed  (black  line,  after  40  hr  low-pass 
filter)  zonal  (upper  panel)  and  latitudinal  (lower  panel)  velocity  components  at  the 
stations  not  used  for  the  reconstruction:  Station-4. 
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Fig.  2c.  OSD  reconstructed  (dashed  line)  and  observed  (black  line,  after  40  hr  low-pass 
filter)  zonal  (upper  panel)  and  latitudinal  (lower  panel)  velocity  components  at  the 
stations  not  used  for  the  reconstruction:  Station- 18. 
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Fig.  2d.  OSD  reconstructed  (dashed  line)  and  observed  (black  line,  after  40  hr  low-pass 
filter)  zonal  (upper  panel)  and  latitudinal  (lower  panel)  velocity  components  at  the 
stations  not  used  for  the  reconstruction:  Station-24. 
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Fig.  3.  TLCS  current  reversals  detected  from  SCULP-1  buoy  trajectories  during  January 
2-7,  1994.  The  black  dots  show  the  starting  positions  of  buoys. 
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98  °W  96°W  94  °W  92°W  90°W 


Fig.  4.  OSD  reconstructed  circulation:  (a)  cyclonic  gyre  on  December  30,  1993,  (b)  total 
reversal  on  January  3,  1994  ,  and  (c)  partial  reversal  on  January  6,  1994. 
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T  (days) 


Fig.  5.  Empirically  (dashed  curve)  and  theoretically  (Poisson  distribution,  solid  curve) 
estimated  recurrence  probabilities:  (a)  Po(T),  (b)  P \(T),  and  (c)  Pi(T)  as  functions  of 
duration  T. 
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Fig.  6.  OSD  reconstructed  mean  circulation  u (x,y)  during  (a)  first  period  (January  21 
-May  21,  1993),  (b)  second  period  (December  19,  1993-  February  12,  1994),  and  (c) 
third  period  (October  5  -  November  29,  1994). 
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Fig.  7.  OSD  reconstructed  first  EOF  mode  ut(x,  v)  during  (a)  first  period  (January  21  - 
May  21,  1993),  (b)  second  period  (December  19,  1993-  February  12,  1994),  and  (c)  third 
period  (October  5  -  November  29,  1994). 
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Fig.  8.  OSD  reconstructed  A\(t)  using  the  LATEX-A  current  meter  mooring  (solid 
curve)  and  the  SCULP-I  drifter  (dashed  curve)  data. 
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Fig.  9.  Time  series  of  (a)  A\{t)  (in:  m/s),  (b)  r/(t) ,  and  (c)  t/ais(0  (in:  m/s).  Eight 
occurrences  of  high  A\ -values  (shaded),  which  represent  the  TLCS  total  current  reversals, 
correspond  to  high  values  of  77  and  Ua is.  Here,  dash-dotted  and  dashed  lines  are  the 
thresholds  for  total  and  partial  current  reversals. 
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Fig.  10.  Morlet  wavelet  spectra  of  (a)  Ua\s(t)  and  (b)  A\(t)  calculated  using  the  data 
collected  during  the  first  period  (January  21  -May  21,  1993).  The  cone  of  influence  for 
A\{t)  is  denoted  by  the  white  curve  above  which  the  edge  effect  (due  to  zero-padding 
ends)  is  greatly  reduced. 


32 


98'W  96°W  94°W  92°W 


Fig.  11.  OSD  reconstructed  circulation  with  northeastly  winds  (around  0.4  m/s)  at  (a) 
surface  and  (b)  15  m  depth  from  combined  SCALP-I  and  Latex-A  observational  data  at 
12:00  GMT  on  February  10,  1994.  Streamfunction  at  15  m  depth  (bold  curves)  shows 
multi-eddy  structure  with  two  cyclonic  eddies  on  the  TLCS.  The  thin  solid  curves  are 
isobaths  at  50,  200,  1000  m  depths. 
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98  °W  96°W  94°W  92°W 


Fig.  12.  OSD  reconstructed  circulation  with  southwestly  winds  (above  3.5  m/s)  at  (a) 
surface  and  (b)  15m  depth  from  combined  SCALP-I  and  Latex-A  observational  data  at 
18:00  GMT  on  January  4,  1994.  Streamfunction  at  15  m  depth  (bold  curves)  shows  an 
anticyclonic  eddy  on  the  TLCS.  The  thin  solid  curves  are  isobaths  at  50,  200,  1000  m 
depths. 
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Fig.  13.  Predicted  (dashed  curve)  and  OSD  reconstructed  (solid  curve)  A\{t)  (in:  m/s) 
during  October-November  1994.  Dash-dotted  and  dotted  lines  indicate  thresholds  for  the 
total  and  partial  TLCS  current  reversals. 
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